Quantum and classical resonant escapes of a strongly-driven Josephson junction 



H. F. Yu, X. B. Zhu, Z. H. Peng, W. H. Cao, D. J. Cui, Ye Tian, 
G. H. Chen, D. N. Zheng, X. N. Jing, Li Lu, and S. P. Zhao 
Beijing National Laboratory for Condensed Matter Physics, 
Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 

Siyuan Han 

Department of Physics and Astronomy, University of Kansas, Lawrence, Kansas 66045, USA 

The properties of phase escape in a dc SQUID at 25 mK, which is well below quantum-to-classical 
crossover temperature T cr , in the presence of strong resonant ac driving have been investigated. The 
SQUID contains two Nb/Al-AlOaj/Nb tunnel junctions with Josephson inductance much larger than 
the loop inductance so it can be viewed as a single junction having adjustable critical current. We 
find that with increasing microwave power W and at certain frequencies v and v/2, the single primary 
peak in the switching current distribution, which is the result of macroscopic quantum tunneling 
of the phase across the junction, first shifts toward lower bias current I and then a resonant peak 
develops. These results are explained by quantum resonant phase escape involving single and two 
photons with microwave-suppressed potential barrier. As W further increases, the primary peak 
gradually disappears and the resonant peak grows into a single one while shifting further to lower 
/. At certain W , a second resonant peak appears, which can locate at very low I depending on 
the value of v. Analysis based on the classical equation of motion shows that such resonant peak 
can arise from the resonant escape of the phase particle with extremely large oscillation amplitude 
resulting from bifurcation of the nonlinear system. Our experimental result and theoretical analysis 
demonstrate that at T <C T cr , escape of the phase particle could be dominated by classical process, 
such as dynamical bifurcation of nonlinear systems under strong ac driving. 
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I. INTRODUCTION 

Devices based on Josephson junctions are not only 
the key elements for the realization of superconduct- 
ing qubitsi but also testbeds for the study of macro- 
scopic quantum phenomena and nonlinear dynamics. In 
these studies, response of the device to resonant or near- 
resonant microwave fields provides insights to under- 
stand their behavior. Numerous works have been re- 
ported in the literature, such as the demonstration of 
energy level quantization in current biased junctions^— 
and Rabi oscillations in phase qubits£~— These experi- 
ments are usually performed at low temperatures (~ 20 
mK) in relatively weak microwave field. Under such con- 
ditions, the junctions behave quantum-mechanically and 
the phase particle may escape from the potential well 
by macroscopic quantum tunneling (MQT) from ground 
state and/or via photon- assisted tunneling (PAT) from 
excited states, depending on the bias current / and the 
microwave frequency v. When the switching current dis- 
tribution P(T) of a current biased junction is measured, 
MQT from ground state manifests as a peak, which we 
call the primary peak, centered slightly below the critical 
current I c of the junction. In the case of the PAT process, 
a resonant peak at a slightly lower current I res could be 
observed in addition to the primary peak in P(I). 

The phase escape of a Josephson junction at low tem- 
peratures exhibits a number of interesting phenomena 
as the microwave power W increases, which can be de- 
scribed as follows. (1) As mentioned above, when the fre- 



quency v matches the bias current dependent energy level 
spacing between the ground and the first excited states 
at I res that is close to I c and W is low, the excited state 
population p\ will be enhanced but remains <C 1. This 
results in a higher escape rate T ~ (1 — p\)To+pxTi > Tq 

at Ires 

since Ti ~ 10 3 ro- Consequently, a resonant peak 

at Ires 

appears in P(I). Because W is low, the system's 
potential and level spacing are essentially unperturbed by 
microwave field, as reflected in the unchanged position of 
the primary peak. (2) When v corresponds to the level 
spacing at a resonance current I res that is well below I c 
and the microwave field is weak, there will be no reso- 
nant peak since Ti is still too low to observe at I res - In 
this case, as W is increased, the position of the primary 
peak will move continuously to lower bias currents due to 
the suppression of the effective potential barrier.— ~— This 
suppression is particularly strong at resonance, leading 
to a significant increase of Ti and thus a visible resonant 
peak. (3) If v matches a fraction of the level spacing and 
satisfies 2nvm — ujq u , where hujQ n is the level spacing be- 
tween the ground state and the nth excited state and m 
is a small positive integer, the multiphoton PAT process 
occurs and the corresponding resonant peak appears. (4) 
When W further increases, the system's nonlinear char- 
acteristics such as nonlinear resonance, bifurcation, and 
chaotic dynamics may emerge. 

The resonant peak has also been observed at higher 
temperatures {e.g. 4.2 K) in the classical regime 
where thermal activation over the top of the barrier 
dominates In this regime, when the driving fre- 
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quency v matches the junction's bias current depen- 
dent plasma frequency uj p (I), the classical resonance 
occurs. The oscillation amplitude of the phase par- 
ticle in the potential well and thus the escape rate 
can be greatly enhanced. Recently, such process and 
their classical interpretation have been studied in ex- 
tended temperature range and have received much at- 
tention due to a series of works by Gronbech- Jensen and 
coworkers These authors demonstrated that previ- 
ously quantum-mechanically explained single- or multi- 
photon resonance, the effective barrier suppression, and 
the Rabi oscillation can also be understood from the clas- 
sical point of view. Interpretations on both the classical 
and quantum bases are also reported^— for the results 
involving multiphoton, multilevel Rabi oscillation, and 
ac Stark shift at high W^— It is found that at high 
W, the predictions of the classical and quantum pictures 
may converged These results indicate that the Joseph- 
son junction system can have a dual character, classical 
and quantum-mechanical, when a microwave field is ap- 
plied. 

In this paper, we investigate the switching current 
distribution and the escape rate of a dc SQUID as a 
function of microwave field strength at low tempera- 
tures. The SQUID contains two Nb/Al-AlO^/Nb junc- 
tions with Josephson inductance much greater than the 
loop inductance so it behaves as a single junction with 
a tunable critical current. We found that as microwave 
power W is increased while frequency v is kept constant, 
the primary peak of P(I) first shifts toward lower bias 
current and then a resonant peak develops. This result 
may find a classical explanation , 1 1 1 1 7 1 1 8 but here we will 
show instead that it is also in accordance with the situa- 
tion (2) discussed above. We further found that the effect 
of two-photon process on escape rate can be observed in 
this situation. Moreover, as W further increases, the 
primary peak will gradually disappear and the resonant 
peak grows into a single one while shifting continuously 
to the lower bias current. At even higher W, a second 
resonant peak appears, which could locate at a very low 
bias current depending on the value of v. Analysis based 
on the classical equation of motion indicates that such 
double-resonant-peak structure with a large separation 
between the two peaks at high W can originate from the 
resonant escape of the phase particle having a very large 
oscillation amplitude resulting from bifurcation of the 
nonlinear system driven by a strong microwave field 



II. EXPERIMENTAL TECHNIQUES AND 
RESULTS WITHOUT MICROWAVE RADIATION 

The sample used in this work was a dc SQUID contain- 
ing two Nb/Al-AlOz/Nb junctions, which was prepared 
at Stony Brook using a self-aligned process.— Each junc- 
tion had a diameter of 1.27 /im (corresponding to an area 
of 1.27 /im 2 ) and the SQUID loop had an area of 220 ^m 2 . 
The device's oscilloscope I-V trace at 25 mK is displayed 
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FIG. 1: (Color online) Switching current distribution P(I) at 
several temperatures between 25 and 800 mK for the sample 
used in this work. The inset shows the oscilloscope I-V image 
at 25 mK. 



in the inset of Fig. 1. Taking I c /2 ~ 1 /LtA, we find the 
Josephson inductance of 330 pH for each junction, which 
is much larger than the SQUID loop inductance of ~ 30 
pH. Hence, the SQUID can be viewed as a single junction 
having a maximum critical current I c . 

According to the RCS J model^S the dynamics of a sin- 
gle junction biased at I can be described by a fictitious 
phase particle with position ip, mass M = C(<I>o/27r) 2 , 
and friction coefficient 1 / RC moving in a washboard po- 
tential 



U(i,tp) = — Ej(iip + cos tp), 



(1) 



in which $o = h/2e is the flux quantum, Ej = I c &o/2ir, 
and i = I / I c . C and R are the capacitance and shunt 
resistance of the junction, respectively. The plasma fre- 
quency of the junction is given by uj p = u>o(l — i 2 ) 1 / 4 , 

where luq = y/2ivI c /C^o- F° r i < 1, the phase parti- 
cle may escape from the potential well cither by ther- 
mal activation (TA) or by macroscopic quantum tunnel- 
ing, resulting in the junction's switching from the zero 
to the finite voltage state. Crossover from MQT to TA 
is known to occur at T cr = %ui p \(l + a 2 ) 1 / 2 — a]/2irkB, 
where a = 1/2Q, and Q = uj p RC is the quality factor. 
The escape rate in the TA regime can be found from 
Kramers' formula^ 



r, 



2^ 



a t exp 



AU \ 



(2) 



where a t = 4/(^/1 + Qk B T/1.8AU + l) 2 is a damping 
dependent factor, and in the MQT regime from^ 



r 9 = 



2tt 



a q exp 



7.2AU 

Tiu)„ 



0.87 s 



(3) 
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FIG. 2: (Color online) Temperature dependence of the mean 
< law > and standard deviation a of the measured P(I) 
(symbols). Solid lines are the predictions from TA and MQT 
theories using the parameters listed in Table I. The experi- 
mental MQT-to-TA crossover temperature T cr — 65 ± 5 mK 
is indicated. 

where a q ~ [1207r(7 .21S.U /hujp)] 1 / 2 . In the above expres- 
sions, the barrier height is given by 

AU(i) = 2Ej[\j\ - i 2 - i arccos(i)] . (4) 

We measured switching current distribution P(I) using 
the time-of-flight technique, as described in our previous 
works . 32 ' 33 The results at some typical temperatures are 
shown in Fig. 1. In our experiment, the bias current 
was ramped up at a constant rate of dI/dt=0A mA/s 
with a repetition rate of 97 Hz. To reduce statistical un- 
certainty each measured P(I) contained 5 x 10 4 escape 
events. The sample was placed in a copper box, which 
was anchored to the mixing chamber (MXC) of an Ox- 
ford Kelvinox MX400 dilution refrigerator. To reduce the 
effects of external noise on switching dynamics to a neg- 
ligible level, a trilayer /li- metal shield was used and the 
circuit was filtered by EMI filters at room temperature, 
cryogenic low-pass filters at 1 K pot, and copper pow- 
der microwave filters at MXC temperature. Since noise 
could cause temperature dependence of the width of P(I) 
to flatten out at T > T cr , which could be mistaken as the 
evidence for MQT, we tested the noise level by applying 
a flux to the SQUID loop to reduce its critical current 
to I c (&) *C / c (0). The result confirmed that noise from 
the environment and measurement circuitry was negligi- 
ble down to 25 mK. (The data presented in this work 
were however obtained with maximum critical current I c 
— I c (0) in the absence of applied magnetic field. The 
plasma frequency at bias current ~ 1.865 jitA, which is 
the position of maximum P(I) at 25 mK in Fig. 1, is lo p 
~ 10 GHz.) 

Fig. 1 shows that as temperature decreases, P{I) nar- 
rows and moves to higher bias current. In Fig. 2, we 
show the temperature dependence of the mean ((I sw )) 
and width (a) of P(I) (symbols), together with the re- 
sults calculated from the TA and MQT theories (lines) 



using the standard conversion between Tt, q and P(/)j24 
From the fit, the sample's parameters I c — 1.957 /iA, 
C = 620 fF, and i? = 300 SI were obtained.^ In Fig. 2, 
we can see that in the temperature range between 1 K and 
T cr ~ 65 mK, a decreases with decreasing temperature 
indicating that TA is the dominant escape mechanism. 
In this temperature range, P(I) depends weakly on C 
and R, so I c can be determined by fitting P(I) using the 
TA theory starting with a reasonable estimate of C and 
R values. The values of C and R were then determined 
using P(I) obtained at T <C T cr because the exponent of 
MQT rate is dependent on C and R. In Table I, we sum- 
marize the sample's key parameters. Taking the P(I) 
peak position at 25 mK in Fig. 1 again, we find that the 
crossover temperature calculated from these parameters 
is T cr = 62.4 mK, which is in good agreement with the 
experimental value indicated in Fig. 2. 



III. RESULTS WITH MICROWAVE 
RADIATION 

Fig. 3(a) shows the switching current distribution P(I) 
at 25 mK in the microwave field with fixed frequency v 
=9.6 GHz and four different power levels. (In this work, 
we always use relative power levels by normalizing to the 
power level at which the resonant peak first appears.) 
Fig. 3(c) shows the corresponding escape rate T. It is 
seen that as W increases, the primary peak of P(I) shifts 
to lower bias current and then a resonant peak develops. 
Further increasing W causes the primary peak to disap- 
pear and resonant peak to grow into a single one (Re- 
sults at higher W are presented in Fig. 4 below). This 
result, though similar to situation (1) discussed in Sec- 
tion I, can be interpreted along a quantum-mechanical 
analysis by Fistul et al£ since the primary peak demon- 
strates a clear leftward shift resulting from effective bar- 
rier suppression^ 

According to the discussion in Ref. |9|, the condition 
{2 / kv/u)q) 5 < Tiuiq/Ej should be satisfied in order for the 
process described in situation (1) to occur. The physics 
behind this is that the microwave-induced excited-level 
population times the escape rate of the level should be 
large. Taking the junction's parameters in Table I, we 
find (27ri//w ) 5 ~ 0.088 and hu /Ej ~ 0.016 so the con- 
dition is clearly not satisfied. On the other hand, we 
calculated the energy levels, and the level escape and 
transition rates of the system using the approach de- 
scribed in Ref. The total escape rate T^o was then 
obtained from the stationary solution of the Master equa- 
tion as first considered by Larkin and Ovchinnikov j 6 ' 37 i 38 



TABLE I: Sample parameters. 
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FIG. 3: (Color online) Switching current distribution P(7) and escape rate T(I) at 25 mK for several microwave powers as 
indicated. Power levels are relative with the curves where resonant peaks start to appear set to dB, and the two dB curves 
in (a) and (b) correspond to -10 dBm and -37 dBm, respectively. Resonant peaks demonstrating resonant quantum phase 
escape can be clearly seen, (a) and (c): v = 9.6 GHz, single-photon process; (b) and (d): v = 4.8 GHz, two-photon process. 
The resonance current I r is indicated by the arrows in (c) and (d). 



We found that for the data without microwave radiation 
(-00 dB) in Figs. 3(a) and (c), there were three levels 
in the potential well and the calculated Tlo agreed well 
with the experimental result. If we assume an unsup- 
pressed potential barrier in the microwave field and use 
Ires =1-79 /i A at the resonant peak, we would have 6 
energy levels in the well. In this case, even assuming 
all population at the first excited level would still lead 
to Tlo about 3 to 4 orders smaller than the measured 
escape rate. From these results, it is clear that a bar- 
rier suppression resulting in a reduction of one or two 
energy levels in the well is necessary to account for the 
experimental data in Fig. 3. 

The calculated level spacing w i/27r between the 
ground and first excited states at I res =1.79 yuA, which 
can also be estimated from ojqi = — 5hui p /36AU), 
is ~ 9.6 GHz, which exactly matches the frequency of 
the applied microwave. This result demonstrates that 
the quantum explanation of the resonant phase escape 
process is applicable. In Figs. 3(b) and (d), we show the 
corresponding results when the microwave frequency is 
reduced to v = (woi/27r)/2 = 4.8 GHz. Similar behav- 
ior can be seen, which is consistent with the quantum 
picture of the resonant phase escape as a result of the 
two-photon process^ 

As W increases, the primary peak in Fig. 3 decreases 
and eventually disappears while the resonant peak grows 
into a single one. The data taken at v = 9.6 GHz are 
plotted in Fig. 4. The most interesting feature of the 



data in this figure is that a second resonant peak devel- 
ops when W is further increased^ Compared with the 
previous case of weak microwave field, the separation be- 
tween the second resonant peak and the first resonant 
peak is much larger than that between the first resonant 
peak and the primary peak. Furthermore, the separation 
is found to depend strongly on the microwave frequency 
v. In Fig. 5, we show the W dependence of P(I) at 25 
mK when v is changed to 7 GHz which is largely de- 
tuned from the level spacing woi/27r ~ 9.6 GHz in the 
vicinity of the bias current where the primary peak is 
located. Therefore, resonant peak originated from PAT 
(e.g. the peaks just below the primary peak in Figs. 3 
and 4) would not occur at low microwave power. As W 
increases, the system gradually enters the strong driving 
regime. The primary peak is observed to shift leftwards 
continuously until a resonant peak at a much lower bias 
current ~ 0.85 /iA develops (note the different horizon- 
tal scales of Figs. 4 and 5). Such distinct results could 
originate from Josephson bifurcation ) 26 ' 27 a phenomenon 
common to strongly-driven nonlinear dynamic systems, 
as discussed below. 



IV. DISCUSSION BASED ON THE CLASSICAL 
NONLINEAR DYNAMICS 

It should be pointed out that the number of levels in 
the potential well at the second resonant peak (at ~ 1.5 
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FIG. 4: (Color online) Switching current distribution P(I) 
at 25 mK with different relative microwave power W. The 
dB curve corresponds to a power level of -10 dBm. The 
microwave frequency is v = 9.6 GHz. It is seen that as W 
increases, the double-peak structure appears and disappears 
twice. The appearance of the two resonant peaks can be ex- 
plained by quantum resonant escape and nonlinear classical 
resonant escape, respectively. 



/iA) in Fig. 4 and the corresponding peak at ~ 0.85 /iA in 
Fig. 5 would be about 20 and over 50, respectively. The 
presence of so many energy levels during resonant escape, 
provided a dramatic barrier suppression did not occur, 
necessitates a many-photon process from the quantum 
mechanical point of view or motion with extremely large 
amplitude oscillations from the perspectives of classical 
dynamics. Since it is not clear whether the treatment of 
Fistul et al. 9 is applicable to such case we provide an ex- 
planation based on the classical description of dynamic 
bifurcation in nonlinear systems for the experimental re- 
sults presented in Figs. 4 and 5. Our analysis shows that 
the experimental results obtained at high W, namely the 
double-peak structure with large peak separation, can be 
caused by escape of the phase particle having very large 
oscillation amplitude resulting from bifurcation of the 
nonlinear system driven by strong microwave field 



A. Nonlinear bifurcation phenomenon: Bistable 
oscillations with different amplitudes 

We start with the normalized equation of motion for 
the phase particle: 



1 



ip- 



Qo 



ip + sin ip = i + i mw cos(rir) 



(5) 



where r = w^b, Qq — u; i?C, i mw — I mw /I c is the nor- 
malized ac current induced by microwave field, and f2 is 
the microwave frequency v normalized to ujq/2tt. At low 
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FIG. 5: (Color online) Switching current distribution P(I) 
at 25 mK with different relative microwave power W. The 
dB curve corresponds to a power level of -15 dBm. The 
microwave frequency is v — 7 GHz. In contrast to the results 
in Fig. 4, the double-peak structure appears and disappears 
only once. Here the resonant peak results from the nonlinear 
classical resonant escape of the phase particle. 



temperatures (T — 25 mK) and strong microwave field, 
the effect of noise current i n due to thermal fluctuations 
on escape rate is negligible and thus omitting i n in Eq. (5) 
is justified. 

Following the analysis of Refs. [l!^ and [H, we use a 
monochromatic ansatz for the solution of the unper- 
turbed Eq. (5): 



with 



ip = tp + 4>{t) 



tpa = arcsin(i) 



(0) 



(7) 



Inserting this ansatz into Eq. (5) and expanding sin</> and 
cos 4> up to the third order yields 



' + i ^ + fll<j) + a<t> 2 + i 



, cos(51t) 



(8) 



in which = (1-i 2 ) 1 / 4 , a = -i/2, and b = -VT~- ^P/Q. 
When the applied microwave frequency fl is close to the 
resonance, we can write f2 = fix + 5£l, where 5il/il <C 1, 
and make the single-mode approximation: 



Asin(f2T + rj) 



(9) 



where rj is a phase constant and A is amplitude of the 
driven oscillation. Inserting this into Eq. (8), one obtains 
the following relation^ 



(10) 
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Solid lines: 9.6 GHz 
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FIG. 6: (Color online) Calculated oscillation amplitude versus 
imw (oc microwave power W) of the phase particle for different 
bias current I and frequency v. The parameters in Table I 
are used. The results show the familiar bifurcation at low I 
and high W , with two oscillation amplitudes A2 > A\ (the 
middle one being unstable for a given curve). 
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FIG. 7: (Color online) Washboard potential U(ip) of Eq. (1) 
with various phase values, remaining barrier height 8U, and 
MQT process indicated. 



B. Resonant phase escape from bifurcated 
oscillation with the larger amplitude: Explanation 
for the results at v = 9.6 GHz 



where x = (36/8fii - 5a 2 /12ftf). 

From Eq. (10), the amplitude A as a function of mi- 
crowave power W oc %^ nw and bias current i can be ob- 
tained since Sfl, x, Qo, and fii only depend on i and/or 
W when microwave frequency VL is constant. Thus, we 
have 



A = A(i,W) . (11) 



In Fig. 6, we present the calculated A(i, W) for three indi- 
cated bias current i at v — 9.6 GHz (solid lines) using the 
parameters listed in Table I. The results show clearly the 
familiar Josephson bifurcation phenomenon! 26 ' 27 Namely, 
as W oc if nw increases, the system can be driven into a 
regime where two dynamic states with different ampli- 
tudes of oscillation A2 > A\ exist, as shown schemati- 
cally in Fig. 7. In Fig. 6, one can see that in the range 
of ij nw considered bifurcation does not occur at bias cur- 
rent I = 1.9 /iA which is slightly below I c (i « 0.971). 
As / decreases, it starts to develop. For lower bias cur- 
rent, threshold for bifurcation shifts to higher W which 
is accompanied by a greater amplitude A2 at the thresh- 
old while A\ has a slow variation (the solid lines in Fig. 
6). The dashed line in Fig. 6 is the result obtained using 
/ = 1.3 /zA and v = 7.0 GHz. In this case, A2 is signif- 
icantly larger when compared to the case of same / but 
higher driving frequency v = 9.6 GHz. This result will be 
shown to be consistent with, and can be used to explain 
the changes from Fig. 4 to Fig. 5 when v decreases. 



We now show that the presence of A2 resulting from 
bifurcation at high microwave power can qualitatively ex- 
plain the double-peak structure with large peak separa- 
tion as presented in Figs. 4 and 5. We first consider the 
four curves with relative power levels of 9, 9.8, 10, and 13 
dB in Fig. 4 with applied microwave frequency v = 9.6 
GHz. In this case, all peaks are located at relatively high 
bias current, so from the results presented in Fig. 6, the 
maximum A2 is seen to be around 1.5, which makes the 
approximation used in deriving Eq. (8) appropriated 

Since our measurement was performed at 25 mK, es- 
cape of the phase particle from the potential via thermal 
activation is negligible. Hence we consider a classical 
resonantly-oscillating particle and its MQT escape from 
the potential well. In general, the particle can escape 
from both A\ and A2 states, thus the MQT rate, in the 
absence of the barrier suppression, can be calculated us- 
ing the WKB method from: 

~ Piv cxp(— 2Si/h) + P2V exp(— 2S2IK) 
~ V 2V exp(-2S 2 /^, (12) 

where px t 2 are the probabilities of the phase particle in 
the A\ and A2 states (pi + P2 = 1) and Si. 2 are the 
corresponding path integral across the remaining tunnel 
barriers. Because S2 -C Si (A 2 ^> Ai) tunneling rate 
from the A\ state is exponentially smaller and the escape 
is dominated by the latter. In an explicit form, we can 
write down the action for the phase particle in the A2 
state as: 
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FIG. 8: (Color online) (a) Calculated 2n5U/hu) p in Eq. (14) 
versus I for the microwave frequency v = 9.6 GHz. Curves 



labeled 1 to 4 have i„ 



2.5x l(T a , 3.7x10" 



4.3x10" 



and 6.5 x 10~ . The parameters in Table I are used in the cal- 
culations. The dotted line is the exponent 7 .2AU (i p ) /Tiuj p {i p ) 
in Eq. (3) where i p is the value corresponding to the primary 
peak without microwave radiation, (b) MQT rate calculated 
from Eq. (12) or (14) using p2 = 10~ 2 and the same v and 
i"rnw parameters from curves 1 to 4 in (a). The line pointed 



by an arrow has different v — 9.4 GHz and i„ 



6.8x 10" 



S-2 



y/U{i,ip) - U(i,ip2a)dif , 



(13) 



in which U is given by Eq. ([I} . The turning point at the 
outer wall of the potential (f2 a = tpo + A2, the point of 
escape ip2b, and the corresponding path of MQT, when 
the phase particle is in the A2 state, are illustrated in 
Fig. 7. 

For the present experimental data, we usually have 
a deep potential well AU [see Eq. (4), for bias current 
well below I c ] and the remaining barrier height SU (see 
Fig. 7) for the particle having A2 oscillation amplitude 
is comparatively small. In this case, it is straightforward 
to show that Eq. (TT21 can be well approximated by 

Fwkb — .P2^exp[-27r<5£//?kjp] , (14) 

where SU is given by 

SU(i, W) = U(i, n-ip )- U(i, ifo + A 2 (i, W)) . (15) 



Various phase locations in Eq. (fT5|) are indicated in Fig. 7. 
In Fig. 8(a), we show the tunneling exponent 2-kSU /Twj p 
of Eq. (U3J) versus I for v = 9.6 GHz as if nw is increased 
from curve 1 to curve 4. Notice that for each curve 
there exists a nearly vertical part, shown as dashed lines, 
which corresponds to the sudden emergence of the A2 
state. This can be understood by considering the case 
for P mw = 3 x 10~ 3 as indicated by the vertical dotted 
line in Fig. 6. At / = 1.3 /iA, only the smaller amplitude 
A\ state exists. The amplitude A\ would increase as I 
increases. At a bias current 1.3 /iA < I < 1.6 /xA, the 
much larger amplitude A2 state emerges which results in 
the sudden decrease of SU(i) and thus a much larger tun- 
neling rate. In Fig. 8(a), it can be seen that the curve 1 
shows a monotonic dependence on /. Starting from the 
curve 2, however, a local minimum develops. This be- 
havior can in principle lead to the double-peak structure 
in P(T) as i'^ nw is increased. 

To compare with T^, which is the MQT rate in the 
absence of microwave radiation, we plot the exponent 
7.2AU(i p )/hujp(i p ) in T q [see Eq. ^ considering Q ~ 

101 as a dotted line in Fig. 8(a), where i p corresponds to 
the position of the primary peak. Notice that the line 
sits well above the wide shoulder part of the solid curves 
resulting from the A2 state of oscillation. Since tunneling 
rate depends exponentially on SU/tiUp our data imply 
P2 ~ exp[2irSU/hoj p — 7 .2AU(i p )/fiojp(ip)] <C 1. Namely, 
the particle spends more time in the Ai state than in the 
A 2 stated 

In Fig. 8(b), we show Twkb calculated from Eq. (fl"2"|) 
or using P2 — 10~ 2 and the same set of parameters as 
those in Fig. 8(a). The choice of P2 around 10~ 2 brings 
the escape rate Twkb for the system in the A 2 state 
of oscillation into the experimentally observable range of 

10 2 to 10 6 sec -1 in Figs. 3(c) and (d), which is deter- 
mined by the bias current sweeping rate of dl / dt = 0.4 
mA/s. It should be noticed that although the calculated 
bias current dependence of Twkb in Fig. 8(b) looks sim- 
ilar to the data presented in Figs. 3(c) and (d) as the 
microwave power increases, the non-monotonic part of 
Twkb(I) covers a much larger range of bias current. 
Furthermore, unlike the data in Figs. 3(c) and (d), the 
hump feature in Fig. 8(b), which lead to the double peak 
P(I) moves much faster to lower currents with increas- 
ing microwave power. These distinct difference led us 
to conclude that the double peak distributions in Fig. 
3 were originated from photon-assisted tunneling rather 
than the dynamic bifurcation under strong ac driving de- 
scribed above. 

However, as the intensity of ac field is increased 
the dominant escape mechanism changes from PAT to 
Josephson bifurcation. Namely, the onset of A2 state re- 
sults in double-peak distributions with large inter-peak 
separation. In Fig. 9, the corresponding switching cur- 
rent distributions P(I) converted from the Twkb{I) pre- 
sented in Fig. 8(b) are plotted. The existence of the 
double-peak structure (curves labeled 2 and 3) can be 
clearly seen. Comparing the results with the experimen- 
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FIG. 9: (Color online) Switching current distributions P(I) 
converted from Vwkb in Fig. 8(b). Curve labels have one- 
to-one correspondence with those in Fig. 8. The data give 
a qualitative explanation for the development of the second 
resonant peak in Fig. 4 (curves with relative powers of 9, 9.8, 
10, and 13 dB). The curve pointed by an arrow is converted 
from that with arrow in Fig. 8(b), and is displayed with its 
magnitude divided by 5. 



tal data in Fig. 4 (namely the curves with relative powers 
of 9, 9.8, 10, and 13 dB), we can see that the shapes of 
P(I) are quite similar and the applied power ratio be- 
tween curve 4 and curve 1 is ~ 2.5, also similar to the 
data displayed in Fig. 4. 



C. Results at v = 7.0 GHz and further discussions 

The results of 9.6 GHz microwave radiation suggest 
that the simple classical model provides a reasonably 
good qualitative explanation for the reappearance of the 
double-peak structure in Fig. 4 as ac driving is intensi- 
fied. Such explanation is expected to work also for the 
data in Fig. 5 with 7-GHz microwave radiation. Unfor- 
tunately, analysis based on Eq. (8) does not work well 
for the v = 1 GHz data since the resonant peak appears 
below 1 fj,A (see Fig. 5). In this case, we can see from 
the dashed line in Fig. 6 that A2 becomes significantly 
greater than tt/2, which makes the expansion of the sin0 
and cos <fi terms in obtaining Eq. (8) inapplicable^ 

Comparing the experimental results in Fig. 4 and 
Fig. 5, there are two clear changes of the bifurcation- 
induced resonant peaks as the microwave frequency v 
decreases. One is the increase of the inter-peak sepa- 
ration of the double-peak distribution and the other is 
the rapid narrowing of the resonant peak width. The 
first change can be easily understood from Fig. 6 that 
for fixed bias current and power, Ai is larger when v is 
lower (compare the solid and dashed lines at I = 1.3 /iA). 
This directly leads to the increase of the inter-peak sep- 



aration for smaller v. It is important to note that such 
trend is basically opposite to that reported in the previ- 
ous work where W cx i 2 mw is smaller and bifurcation is 
not involvedi 17 ' 18 

The rapid narrowing of the peak can also be explained 
by our model. To demonstrate this, we plot in Fig. 8(b) 
the result with v — 9.4 GHz (the top curve indicated by 
the arrow), a value smaller than 9.6 GHz used in the fig- 
ure but not far away so that the approximation used for 
Eq. (8) is still valid. We have raised i^ nw cx W slightly so 
that the vertical part of the curve sits at the same bias 
current as that of curve 4. We can see that the calculated 
^wkb increases quickly, roughly by one order of magni- 
tude when v is decreased only by 0.2 GHz. Such increase 
would influence the peak width significantly since, unlike 
the usual case such as those depicted in Figs. 3(c) and (d) 
where escape rate changes continuously, the bifurcation 
related escape rate versus dc bias current makes a sud- 
den jump (dashed parts of the curves in the figure) when 
Ai appears. While the rate of escape from the smaller 
A\ amplitude oscillation is negligible, it increases rapidly 
when the A2 oscillation appears, which results in switch- 
ing to occur in a very small range of bias current and 
therefore a much narrower resonant peak. Apparently, 
the lower the v is, the higher level the escape rate will 
suddenly jump to, and the narrower the resonant peak 
will be. The curve pointed by the arrow in Fig. 9 is the 
resulting distribution that appears considerably sharper 
than curve 4 in the same figure. 

It should be pointed out that from our simple model 
the "primary" and "resonant" peaks in Fig. 9 have differ- 
ent meanings as conventionally understood from the data 
such as those in Fig. 3. According to the discussion on the 
results presented in Fig. 8, all peaks in Fig. 9 are caused 
by tunneling from the A2 state of the resonantly driven 
phase particle. The double-peak structure is a result of 
combined effect of non-trivial bias current dependence of 
the potential well depth U(i,ir — ipo) — U(i,ipo) and the 
amplitude A 2 (i, W) at a given W. According to Eqs. (1) 
and (11), both quantities decrease as i increases but at 
different rates. This gives rise to the local minima and 
humps in Fig. 8, and the "resonant" peak in Fig. 9. 

We note that quantitatively the experimental data and 
calculation based on Josephson bifurcation do not agree 
well, thus further improvements of the model are neces- 
sary. Among them, including higher order terms in the 
expansion of sin (f> and cos <fi should be considered for in 
order to describe the data taken at v = 7.0 GHz4i We 
emphasize that Josephson bifurcation is used to explain 
the appearance of the second resonant peak in Fig. 4 
for data taken at 9.8 dB and 10 dB. The double-peak 
structure at much weaker ac driving, from about dB 
to 3 dB, is due to the quantum resonant phase escape 
process accompanied by barrier suppression;*^^ which is 
not considered in the present model. 
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V. SUMMARY 

We presented a systematic study of resonant phase es- 
cape in tilted washboard potential of a Josephson junc- 
tion (a dc SQUID with small loop inductance) at 25 mK 
in the presence of strong ac driving (microwave radia- 
tion). The device was well-characterized using the TA 
and MQT theories in the absence of the microwave field. 
For weak ac driving i mw <C 1, we observed results that 
are interpreted by quantum resonant escape of the phase 
involving single- and two-photon absorptions with a sup- 
pressed potential barrier. At larger i mw , a second res- 
onant peak, which is well separated from the first one, 
would appear. The peak could locate at very low bias 
current depending on the power and frequency of the mi- 
crowave applied. We proposed a model based on the clas- 
sical equation of motion to interpret the data. Our results 
indicated that at large i mw , the phase particle enters a 
bistable state due to bifurcation of the nonlinear system 
and the oscillation state with larger amplitude leads to 
the resonant peaks in the switching current distributions 
locating far below the system's critical current. These 
results are useful for the further studies of the nonlin- 



ear response of Josephson junctions to strong ac driving 
and for the detection of the qubit quantum states us- 
ing a dc SQUID. In addition, the results and analysis 
showed that for junctions with strong ac driving, even 
at T <C T cr , the quantum-to-classical crossover tempera- 
ture, escape of the phase particle could be dominated by 
classical processes. Our work therefore provides further 
evidence of the dual character, classical and quantum me- 
chanical, of the Josephson junction system at T <C T cr 
in an ac driving field. 
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